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Abstract 

In this article, we calculate the strong coupling constant qantt and study 
the strong decay A"'"^ — s- pn with the light-cone QCD sum rules. The nu- 
merical value of the strong coupling constant g/\Nn is consistent with the 
experimental data. The small discrepancy maybe due to failure to take into 
account the perturbative 0{as) corrections. 

PACS numbers: 13.30.-a; 13.75. Gx 

1 Introduction 

The A(1232) resonance dominates many nuclear phenomena at energies above the 
pion-production threshold and plays an important role in the physics of the strong 
interaction. It is almost an ideal elastic nN resonance, and decays into the nucleon 
and pion (A Nn) with the branching fraction about 99%. The only other (elec- 
tromagnetic) decay channel (A — > A^7) contributes less than 1% to the total decay 
width [1]. There is a very small mass gap (less than 300MeV) between the A and 
the nucleon, and the A(1232) is taken as an explicit dynamical degree of freedom in 
the heavy baryon chiral perturbation theory pj. 

In this article, we calculate the strong coupling constant gANn with the hght-cone 
QCD sum rules, and study the decay width PA^iVTr- The strong coupling constants 
of the octet baryons with the vector and pseudoscalar mesons Qnnv and Qmnp have 
been calculated with the light-cone QCD sum rules [3]. The light-cone QCD sum 
rules carry out the operator product expansion near the light-cone ~ instead 
of the short distance x ~ while the nonperturbative hadronic matrix elements 
are parameterized by the light-cone distribution amplitudes instead of the vacuum 
condensates [H [5] . The nonperturbative parameters in the light-cone distribution 
amplitudes are calculated with the conventional QCD sum rules and the values are 
universal |6j. 

The article is arranged as: in Section 2, we derive the strong coupling constant 
gANn with the light-cone QCD sum rules; in Section 3, the numerical result and 
discussion; and Section 4 is reserved for conclusion. 
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2 Strong coupling constant ^aa^tt with light-cone 
QCD sum rules 

In the following, we write down the two-point correlation function n^(p, g), 

Il^{p,q) = ijd^xe-''^-^{0\T{jp{0)J^{x)}\7T{p)), (1) 
Jp{x) = e''''%^(x)C7^Mb(x)757''4(a;) , 

J,ix) = e''V^{x)C-f^u,{x)u,ix) , (2) 

where the baryon currents Jp{x) and Jfj.{x) interpolate the octet baryon p and de- 
cuplet baryon respectively jY] , the external state vr has the four momentum 
with = . The general form of the proton current can be written as |8] 

Jp{x,t) = Eabc { [ul{x)Cdb{x)\ j5Uc{x) + t [u^ {x)C'j5db{x)] Uc{x)} , 

in the limit t = —1, we recover the loffe current. If we retain the additional pa- 
rameter t and choose the ideal value, the sum rule maybe improved, in this article, 
we choose the loffe current for simplicity. The correlation function n^(p, g) can be 
decomposed as 

Ilsig^ + 11^2 + nij3 pfq^ + nB4cr„/3p"g^g^ + 

nciTM + qhti + ncs i/it, + ^ci(^^Jiuo,(}Yl5p'^q'^ (3) 

due to the Lorentz invariance, where the 11 and Ilj are Lorentz invariant functions 
of p and q. In this article, we choose the tensor structure (TajsP^^q^P^ for analysis. 

The strong coupling among the A, p and tt can be described by the following 
chiral Lagrangian [2], 

C{x) = g^M.[^^{x)d^'K{x)N{x)+N{x)d^'K{x)^^'{x)\ . (4) 

Basing on the quark-hadron duality [6], we can insert a complete series of inter- 
mediate states with the same quantum numbers as the current operators Jp{x) and 
J^(a;) into the correlation function n^(p, g) to obtain the hadronic representation. 
After isolating the ground state contributions from the pole terms of the baryons p 
and A, we get the following resullH , 

^ ^ (0|Jp(0)|iV(g + p))(iV(g + p)|A(g)7r(p))(A(g)|J^(0)|0) 
'^^'^^ {M^,-{q + P?]{Ml-q^} 

= {M^p - (g + {Ml - g^} i^""^^"^^^'^ + ■■■) + ■■■, (5) 



^In the first version of tliis article(lar Xiv:0707.3736 |, the numerical factor -i is missed in 



3 -CaisPaQfBPfj, and we obtain too smaU value for the strong coupling constant qan-k to accom- 
modate the experimental data. 
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where the following definitions have been used, 

(0|J,(0)|iV(p)) = xMp,s), 

(0|J^(0)|A(p)) = XaU,{p,s), 

s 

{N(g')\A(g)n{p)) = ig„NjJ(q' ,s')U„(<i,s)jf , (6) 

the last identity corresponds to the phenomenological Lagrangian in Eq.(4). 

The current J^{x) couples not only to the isospin / = | and spin-parity = \ 
states, but also to the isospin / = | and spin-parity = \ states. For a generic 
i resonance A* [9], 

(0| J,(0)|A*(p)) = A.(7, - 4^)f/*(P, s) , (7) 

where A* is the pole residue and M^: is the mass. The spinor U*{p,s) satisfies 
the usual Dirac equation — M^)U*{p) = 0. If we take the phenomenological 
Lagrangian, 

C(x) = gj^,N^{A*{x)N{x)7T{x) + N{x)A*{x)n{x)} , (8) 

which corresponds to {N{q')\A*{q)'K{p)) = gA*NTrU{q' , s')U*{q, s), the contributions 
from the | states can be written as 

+ ■■■ 

= IlDqf)f, + IlEiqf^ + IIe2 qil^i + ^es pflf, + ^EACfapp'^q'^q^i + 

lipiiiJi + qhil + n^s p^ifi + ^Fi^iiuapi^'ibP'^q^ , (9) 

where the Ilj are Lorentz invariant functions of p and q. If we choose the tensor 
structure (Ta/sPaQ/sPfi, the A* has no contaminations. 

In the following, we briefiy outline the operator product expansion for the corre- 
lation function 11^ (p, q) in perturbative QCD theory. The calculations are performed 
at the large space-like momentum regions {q + p)"^ ^ and <C 0, which corre- 
spond to the small light-cone distance ~ required by the validity of the operator 
product expansion approach. We write down the "full" propagator of a massive light 
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quark in the presence of the quark and gluon condensates firstly [U [6p, 



^ / \ i^ab ^ab^q ^ab / _ \ , ^^ab / _ \ ^ab'^ / _ ^ \ 

Dabix) = — TT-r ?r-?r (w) H ^aKQQ) KQQs'y^Q) 

^ ' Atx^x^ 12 48 192 ^' 



idabX^ 

+ -^^mg{qgsaGq) ^ 



1 

dv [(1 - v)gsG^^{vx) ^af"" + vg^G ^^{vx)a^'' 



IGtt^x^ 

+ ■■■, (10) 

then contract the quark fields in the correlation function n^(p, g) with the Wick 
theorem, and obtain the following result, 

n^.(p,g) = '2ieabc^a'b'c> j d^xe'"^"" 

{Tr [^^Shb'i-x)j,GSl,{-x)C] 757"(0|4(0)«e'(x)|7r(p)) 

-2j,r{o\dMM^)Hp)KcsU-^)CiaSbA-x)} . (11) 

Perform the following Fierz re-ordering to extract the contributions from the two- 
particle and three-particle yr-meson light-cone distribution amplitudes respectively, 

1 1 

<la{0)qp{x) = -—Sab5a(iqix)q{0) - —6ab{7'')a/3q{x)-f^q{0) 

-^^ab{(r^'')apq{x)(r^^q{0) 

+ J^^ab{l''l5)af3q{x)jf,-f5q{0) 

+Y^Sab{ij5)apq{x)i-f5q{0) , (12) 

q:{0)(f^{x)Gt{vx) = -^6^pqix)Gxr{vx)q{0) - ^{Y)apq{xh,Gxr{vx)q{0) 

-^{cr^'')af^q{x)crf,^GxT{vx)q{0) 
o 

+^h^'y5)apq{xhf,-f5Gxr{vx)q{o) 

+ ^{i'l5)a/3q{x)i-f5Gxr{vx)q{0) , (13) 

and substitute the hadronic matrix elements (such as the (0|tl(x)7^75(i(0)|7r(p)), 
{0\u{x)gsaf,y^5Gaf3{vx)d{0)\7i{p)), (0|M(x)cr^,,75(i(0)|7r(p)), etc.) with the correspond- 
ing TT-meson light-cone distribution amplitude^, finally we obtain the spectral den- 
sity at the coordinate space. Once the spectral density in the coordinate space 



■^One can consult the first article of Ref. [4] and the second article of Ref.[6j for the technical 
details in deriving the full propagator. 



^ In calculations, we have used the relations cr^^ = —■k^^i,uaf5cr°'^lb and G^^ = ^e^^i/a/jG"^. 
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is obtained, we can translate it into the momentum space with the D = A + 2e 
dimensional Fourier transform, 



4n = — / duu(p^{u) — -- / duuA{u^ 



3vr2 7o '{-Q'Y 27r2 J, ^ ' {-Q 

U asGG [\ , r(2) 



2U 



36 W ' 7o ' ' {-Q'Y 

1 rl rl-aa 



dvv / dag / da^ 
io Jo 



r(2) 



2^2 



-Q 

f^ml /"^""s r(i 



Jo Jo Jo \^ ) 

[{1 - 9v)V^_ - 4(1 - 2v)A\i - 4(1 - v)A^] (a„ a^, 1 - a„ - 

■^'^ ^ ^ dv f dan f da„ f da- ^ ^ ' 



+ F± + (1 - 2v){A\\ + A±)] {a, ag,l-a- ag) 

rft;(l -V) \ dag \ dp "'^ /_n2u l«=i-(i-f)"s 



Jo -/o A Jo ( 

[V^l + + (1 - 2t;)(A|j + A^)] (a, 1 - a - /?) , (14) 

where Q^i = q^i + upfj, and = (1 — u)q'^ + u{p + g)^ — u{l — u)rn^. The e is a 
small positive quantity, after taking the double Borel transform, we can take the 
limit e — > 0. 

There is no contribution from terms of the form {qq)(f)n{u), while there is rather 
large contribution from that terms in the sum rules for the strong coupling constant 
gNN-K, see the article "V. M. Braun and I. E. Filyanov, Z. Phys. C44 (1989) 157" 
in Ref. [1]. If we replace the decuplet baryon current J^{x) with the octet baryon 
current Jn(x) (interpolating the neutron) and study the strong coupling constant 
dNN-K, the Feynman diagrams are quite different. Our mathematica code can be used 
to calculate the strong coupling constant gNNn and produce the terms {qq)4'-Kiu). 

The decuplet baryon current J^{x) and octet baryon current Jp{x) have the 
Dirac structures ® 1 and 7^ ® 7^75 respectively, where ® stands for the u quark 
fields. The Dirac structure 7^75 corresponds to the twist-2 light-cone distribution 
amplitude (pniu). If we replace one of the "full" u quark propagators with the 
quark condensate (qq), the terms {qq)4>n{u) in the correlation function 11^ (p, q) have 
the Dirac structures j^jalx or 'jfj.'yx'ja, which are chiral even, because only the 
perturbative part of the other "full" u quark propagator has contribution. It is not 
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unexpected that they have no contribution to the chiral odd structure cfaiBPaQ/sPfi- 

The hght-cone distribution amphtudes (/>7r(ii), A{u), (f>3Tr{ai), A^{ai), A\\{ai), 
V±{ai) and V||(aj) of the vr meson are presented in the appendix [TU], the nonper- 
turbative parameters in the hght-cone distribution amphtudes are scale dependent, 
in this article, the energy scale is taken to be /i = 1 GeV. The contributions propor- 
tional to the G^i, can give rise to three-particle (and four-particle) meson distribution 
amplitudes with a gluon (and quark- ant iquark pair) in addition to the two valence 
quarks, their corrections are usually not expected to play any significant role^. In 
this article, we take them into account for completeness. 

Taking double Borel transform with respect to the variables Qi = —q^ and Q\ = 

-ip + qf respectively (i.e. ^^^,^^^^,J^t)Ql+uQl\^ ^^^"'^^^(^ - ^o), 

= M'^+th ^0 ~ M^V^ ) ' then subtract the contributions from the high 
resonances and continuum states by introducing the threshold parameter sq (i.e. 
M^" — i> Yjn] /o*° dss"'~^e~^), finally we obtain the sum rule for the strong coupling 

^For examples, in the decay B — s- XcqK, the factorizable contribution is zero and the nonfac- 
torizable contributions from the soft hadronic matrix elements are too small to accommodate the 
experimental data [llj; the net contributions from the three- valence particle light-cone distribu- 
tion amplitudes to the strong coupling constant gD^iD'K are rather small, about 20% [12]. In 
Ref. [13], we observe that the contributions from the three-particle (quark-antiquark-gluon) light- 
cone distribution amplitudes are less than 5% for the strong coupling constants go'DT- In this 
article, the contributions from the three-particle light-cone distribution amplitudes are about 10%. 
The contributions from the three-particle (quark-antiquark-gluon) distribution amplitudes of the 
mesons are always of minor importance comparing with the two-particle (quark-antiquark) distri- 
bution amplitudes in the light-cone QCD sum rules. In our previous work, we also study the four 
form- factors /i(Q^), /2(Q^), gi{Q^) and g2{Q^) of the S ^ rt with the light-cone QCD sum rules 
up to twist-6 three-quark light-cone distribution amplitudes and obtain satisfactory results |14j . 
In a word, we can neglect the contributions from the valence gluons and make relatively rough 
estimations in the light-cone QCD sum rules. 



6 



constant (^aatti 



XpXA 



exp 



Ml 
Ml M| 



Ml uq{1 - uo)ml\ f Mo 



M2 



87r2 



36 TT 



144' TT 



M2 



uo—au 

UQ 



87r2 



4(i-2 ^°:"" Mji+4(i- ^°:"" )A, 



an 



an 



(q^u) ctg) 1 a^i ag^ 



;M^Eo{x)Uml 



l-UO 



1-MO 



da,, 



UQ—ag 



dan 



da 



MO 

1 

a„ 



da 



V\\ + V^ + {l-2'^^-^){A\\+A, 



an 



+—M'Eo{x)Umi{l-uo) dag— d(3 



TT 



l-UQ 



1-/3 



V|| + + (1 - 2 



an 



, (15) 



where 



En{x) 



X 



l-(l + x + - + --- + -)e-% 
M2 ■ 



In the following, we present an ansatz for the spectral density at the level of 
quark-gluon degrees of freedom [15], [16]. Firstly, we perform a double Borel trans- 
form for the correlation function (which is denoted as du-rrt — 2[^}J'^'^\r^i , ^aio 
symbolically) with respect to the variables Q\ and Q\ respectively, and obtain the 
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result, 



M2{2-a) 

Ml Ml 



du- 



T{a)f{u) 



exp 



u{l - u)ml + (1 - u)Ql + iiQ2]. 



M2 



(16) 



where the f{u) stand for the hght-cone distribution amphtudes, a < 2, uq = j^ij^j^i , 



2 M?M2 



1 ^"2 

M'(+Ml ■ 



Then we introduce the corresponding spectral densities p(si, S2), 



M2(2-a) exp 



2 n 



%(1 - 'Wo)"^^ 



M2 



f{uo) 

S2 



Ml 



Ml 



ds2 exp 

and take a replacement Ml M| — » 

cis2 exp {-siCTi - 520-2} p(si, S2) 
exp {-^0(1 - UQ)ml{ai + aa)} 



(17) 



/(«o) 



(0^1 + 0^2 

/(^o) 

r(2-a; 



\2-a 



POO 

/ (iAA^-"exp {- [mo(1 - uo)ml + X] (ai + ^2)} . (18) 
Jo 



Finally we take a double Borel transform with respect to the variables ai and (T2 
respectively, the resulting QCD spectral densities read 



r(2 - a] 



{si - uo{l - uo)miy " 6{si - S2 



(19) 



The threshold parameter sq is taken as Sq = max(s5,S2), where the and 
are the threshold parameters for the channels 1 and 2 respectively. The quantity 
uq{1 — UQ)m^ is tiny and can be safely neglected. Our approach (i.e. performing a 
double Borel transform and taking a replacement M^" — > f^^ dss"'~^e~^ .) is 
an indirect way to obtain the same results. 



3 Numerical result and discussion 

The input parameters are taken as niu = nid = (0.0056 ± 0.0016) GeV, /tt = 
0.130 GeV, = 0.138 GeV, A3 = 0.0, /g^ = (0.45 ± 0.15) x lO^^GeV^ us = 
-1.5 ± 0.7, = 0.2 ± 0.1, 02 = 0.25 ± 0.15, ai = 0.0, = 10.0 ±3.0 [H [TDl [H], 
(qq) = -(0.24 ± 0.01 GeV)3, {^) = (0.33 GeV)^ [6\, Mp = 0.938 GeV, Ma = 
1.232 GeV, Ap = (2.4 ± 0.2) x 10 GeV^ and Aa = (3.0 ± 0.2) x 10"^ GeV^ ^7]. 
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In this article, we neglect the perturbative 0{as) corrections to the strong cou- 
pling constant gANw, and take the values of the pole residues Xp and Aa without 
perturbative 0{as) corrections for consistency. 

In calculation, we observe the main uncertainties come from the two parameters 
02 and in the two-particle light-cone distribution amplitudes, as the dominant 
contributions come from the two-particle light-cone distribution amplitudes (pniu) 
and A{u), the contributions from the terms involving the three-particle (quark- 
antiquark-gluon) light-cone distribution amplitudes are of minor importance, about 
7% of the contribution from the term |^M^i?i(x)/7r07r(^^o)- The uncertainty of 
the parameter 02 obtained in Ref.flO] is very large, in this article, we take smaller 
uncertainty, say 30% (i.e. 02 = 0.25 ± 0.08), which is the typical uncertainty in the 
QCD sum rules. 

The values of the vacuum condensates have been updated with the experimental 
data for the r decays, the QCD sum rules for the baryon masses and analysis of 
the charmonium spectrum [HI [191 [2Q], in this article, we choose the standard (or 
old) values to keep in consistent with the sum rules used in determining the non- 
perturbative parameters in the light-cone distribution amplitudes. 

The threshold parameter sq is chosen to be sq = (3.4±0.1) GeV^ to avoid possible 
contamination from the contribution of the P33 baryon A(1920) in the pn'^ scattering 
amplitude |2l]. Furthermore, it is large enough to take into account the contribution 
of the A(1232). However, the interpolating current J^(x) has nonvanishing coupling 
with the isospin ^ = ^ and spin J = ^ states, the contribution from the 5*31 state 
A(1620) is included in if the A(1620) has negative parity [T]. We choose the tensor 
structure (Ta/sPalpPiJ, to avoid the contamination. 

The Borel parameters are chosen as ^ = ^ and = ]^i^]^2 = (2.0 — 

3.0) GeV^, in those regions, the value of the strong coupling constant gANn is rather 
stable with the variation of the Borel parameter M^, which are shown in Fig.l. 

Taking into account all the uncertainties, finally we obtain the numerical results 
for the strong coupling constant gANn, which are shown in Fig.l, 



gAN^ = (13.5 ± 7.2) GeV"^ 
gANn = (13.5 ± 5.4) GeV"^ 



(20) 



for the parameter 02 = 0.25 ± 0.15 and 02 = 0.25 ± 0.08 respectively. 

The strong coupling constant g^Nn has the following relation with the decay 
width Fa^attt, 



A^Ntt 



ss' 



^[Mj - (M, + m^nM'i - (M, , 

Pcm - ■ ^ ' 

If we take the experimental data as input parameter, Fa-^attt =118 GeV pj, we can 
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M'[GeV'] M^[GeV^] 



Figure 1: Qantt with the Borel parameter M^, A for 02 = 0.25 ± 0.15 and B for 
as = 0.25 ±0.08. 

obtain the value gANw ~ 15.6 GeV~^, our numerical result g^Nn = (13.5±5.4) GeV^^ 
is rather good. 

In the region = (2.0 - 3.0) GeV^ ~ 0.10 - 0.13 [20]. If the radiative 

0{as) corrections to the leading perturbative terms are companied with large nu- 
merical factors, just like in the case of the QCD sum rules for the mass of the proton 

1 + + 1e)^^ - 1 + (0.54 - 0.65) , (22) 

the contributions of the order 0{as) are large, neglecting them can impair the pre- 
dictive ability. Furthermore, the pole residues Xp and Aa also receive contributions 
from the perturbative 0{as) corrections, if we take them into account properly, we 
can improve the value of the strong coupling constant Qan-k- 



4 Conclusion 

In this article, we calculate the strong coupling constant qan-k and study the strong 
decay A"*""^ — > pn^ with the light-cone QCD sum rules. The numerical value of 
the strong coupling constant qan-k is consistent with the experimental data. The 
small discrepancy maybe due to failure to take into account the perturbative 0{as) 
corrections. 
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Appendix 

The light-cone distribution amphtudes of the tt meson are defined by fLO 

1 C 2 2 

' mix 







+ y^ml^ / due-'-^-^B{u) 
2 p- X Jq 

,2 n 



{0\u{x)t-f,d{0)\7r{p)) = JlII^ j due-'^P-^<j),{u) , 

rriu + md Jo 



.2 /-l 



f m 

(0|M(x)a^^75c/(0)|7r(p)) = i{Xf,p^ - x^p^) "" J I due-'''P-'''(f)^{u) 

6(m„ + md) Jo 

{0\u{x)af,^-f59sGc,f3{vx)d{0)\7T{p)) = /s^ {{p^p^g^f^ - PuPaQ^p) - {Pf^PisgL 
{0\u{x)jf,j5gsGa(3{vx)d{0)\7r{p)) = /^m^p^?"^^ 



p ■ X 

+f^ml{p^g^^-Po,g^^) 
{0\u{x)'yi,igsGaf^{vx)d{0)\7T{p)) = /^m^p^^"^^ 



p ■ X 



ai e 



-ip-x(au+vag) 



+f^m^{pf3g^^~Pagf3^) 

j Paiyx(ai)e-'P "(""+''"«) , (23) 



where g^^^ = gf,y - p^^-^+p-^^ ^ (5^^ = ie^^.a/jG"'^ and Vat = daudaddag6{l - Uu - 

OLd - OLg)- 
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The light-cone distribution amphtudes of the tt meson are parameterized as [TO] 



A{u) 



9{u 



u 



(j)j,{u) = 6u{l — 



u = 



10' 
1 



+ |-3?73CJ3 

6n(l - n) |l - 

+ A3(a„ - ad) + uj3^{7ag - 3)| , 
= UOauadttg (voo + Vio{3ag - 1)) , 



3 
5' 



Q(0 



037r(ai) = 360auadal \ 1 
A\\{ai] 



= 120auC(dGtgO'lo{c(d — « 

= -30a^ {/ioo(l - ctg) + V K(l - ttg) - 6a„ad] 



■10 



a 



30a^(a„ - ttd) 1^00 + hiag + ^^io(5ag - 3) 
f 16 24 20 

QU{1 - ^) ^ + + 20r/3 + y774 

1 T '7 1 M -3 



+ 



+ 



i i 7 10 

~ 15 ^ 16 ~ 27 ^^'^^ ~ 27^^ 
11 4 1 a 



7^2 — 



135 



Q(0 



-^02 + 21r]iUi 



^ 210 J J k " ; 

{2ii3(10 - 15m + 6^2) logM + 2^3(10 - 15ti + Qu^) \ogu 
+uu{2 + 13mm)} , 



B{u) = g{u) - (l)^{u) 



(24) 
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where 



h - _ ^4 
"-oo — ''^oo — ~y 1 

21 9 
21 

^10 = -^fli^i , 
o 

h ^ 3 

7 3 

^10 = -^Vi^i + , 

18 20 
5^2 = 1 + y a2 + 6O773 + y ?74 , 

9 

94 = -^«2-6r73aj3, (25) 



^ = 2^ — 1, and 6*2^(0; C*! (0; C2^(^) are Gegenbauer polynomials, 773 



Acknowledgments 

This work is supported by National Natural Science Foundation, Grant Number 
10405009, 10775051, and Program for New Century Excellent Talents in University, 
Grant Number NCET-07-0282. The author would like to thank Prof.T. Huang and 
Prof.V.M.Braun for valuable discussions and comments. 



References 

[1] W.-M. Yao et al, J. Phys. G33 (2006) 1. 

[2] T. R. Hemmert, B. R. Holstein and J. Kambor, J. Phys. G24 (1998) 1831. 

[3] T. M. Aliev, A. Ozpineci and M. Savci, Phys. Rev. D64 (2001) 034001; T. Doi, 
Y. Kondo and M. Oka, Phys. Rept. 398 (2004) 253; T. M. Aliev, A. Ozpineci, 
S. B. Yakovlev and V. Zamiralov, Phys. Rev. D74 (2006) 116001; Z. G. Wang, 
Phys. Rev. D75 (2007) 054020; and references therein. 

[4] I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, Nucl. Phys. B312 (1989) 
509; V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B345 (1990) 137; V. L. 
Chernyak and A. R. Zhitnitsky, Phys. Rept. 112 (1984) 173; V. M. Braun and 
I. E. Filyanov, Z. Phys. C44 (1989) 157; V. M. Braun and I. E. Filyanov, Z. 
Phys. C48 (1990) 239. 



13 



V. M. Braun, hep-ph/9801222 P. Colangelo and A. Khodjamirian, 
hep-ph/0010175] 

M. A. Shifman, A. 1. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 
385, 448; L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 
1; S. Narison, QCD Spectral Sum Rules, World Scientific Lecture Notes in 
Physics 26 (1989) 1. 

B. L. loffe, Nucl. Phys. B188 (1981) 317; B. L. loffe and A. V. Smilga, Nucl. 
Phys. B232 (1984) 109; V. M. Belyaev and B. L. loffe, Sov. Phys. JETP 56 
(1982) 493; V. M. Belyaev and B. L. loffe, Sov. Phys. JETP 57 (1983) 716. 

V. Chung, H. G. Dosch, M. Kremer and D. SchoU, Nucl. Phys. B197 (1982) 
55; H. G. Dosch, M. Jamin and S. Narison, Phys. Lett. B220 (1989) 251. 

V. M. Braun, A. Lenz, G. Peters and A. V. Radyushkin, Phys. Rev. D73 (2006) 
034020. 

P. Ball, JHEP 9901 (1999) 010; P. Ball and R. Zwicky, Phys. Lett. B633 (2006) 
289; P. Ball and R. Zwicky, JHEP 0602 (2006) 034; P. Ball, V. M. Braun and 
A. Lenz, JHEP 0605 (2006) 004. 

L. Li, Z. G. Wang and T. Huang, Phys. Rev. D70 (2004) 074006; B. Mehc, 
Phys. Lett. B591 (2004) 91. 

Z. G. Wang, J. Phys. G34 (2007) 753. 

Z. G. Wang, Nucl. Phys. A796 (2007) 61. 

Z. G. Wang, J. Phys. G34 (2007) 493. 

V. A. Beilin and A. V. Radyushkin, Nucl. Phys. B260 (1985) 61; P. Ball 
and V. M. Braun, Phys. Rev. D49 (1994) 2472; V. A. Nesterenko and A. V. 
Radyushkin, Phys. Lett. B115 (1982) 410; H. Kim, S. H. Lee and M. Oka, 
Prog. Theor. Phys. 109 (2003) 371. 

Z. H. Li, T. Huang, J. Z. Sun and Z. H. Dai, Phys. Rev. D65 (2002) 076005. 

V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Riickl, Phys. Rev. D51 
(1995) 6177. 

[18] B. L. loffe and K. N. Zyablyuk, Nucl. Phys. A687 (2001) 437; B. V. Geshken- 
bein, B. L. loffe and K. N. Zyablyuk, Phys. Rev. D64 (2001) 093009; B. L. 
loffe and K. N. Zyablyuk, Eur. Phys. J. C27 (2003) 229; K. Zyablyuk, JHEP 
0301 (2003) 081. 

[19] B. L. loffe. Prog. Part. Nucl. Phys. 56 (2006) 232; and references therein. 



14 



[20] M. Davier, A. Hocker and Z. Zhang, Rev. Mod. Phys. 78 (2006) 1043; and 
references therein. 

[21] R. A. Arndt, W. J. Briscoe, I. I. Strakovsky and R. L. Workman, Phys. Rev. 
C74 (2006) 045205. 



15 



